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ABSTRACT 

The Infrared Spectrograph (IRS) on board the Spitzer Space Telescope observed about 15 000 objects 
during the cryogenic mission lifetime. Observations provided low-resolution (R = A/AA « 60 - 127) 
spectra over « 5 - 38 pm and high-resolution (R « 600) spectra over 10 - 37 pm. The Cornell Atlas 
of Spitzer/IRS Sources (CASSIS) was created to provide publishable quality spectra to the community. 
Low-resolution spectra have been available in CASSIS since 2011, and we present here the addition of the 
high-resolution spectra. The high-resolution observations represent approximately one third of all staring 
observations performed with the IRS instrument. While low-resolution observations are adapted to faint 
objects and/or broad spectral features (e.g., dust continuum, molecular bands), high-resolution observa¬ 
tions allow more accurate measurements of narrow features (e.g., ionic emission lines) as well as a better 
sampling of the spectral profile of various features. Given the narrow aperture of the two high-resolution 
modules, cosmic ray hits and spurious features usually plague the spectra. Our pipeline is designed to 
minimize these effects through various improvements. A super sampled point-spread function was cre¬ 
ated in order to enable the optimal extraction in addition to the full aperture extraction. The pipeline 
selects the best extraction method based on the spatial extent of the object. For unresolved sources, the 
optimal extraction provides a significant improvement in signal-to-noise ratio over a full aperture extrac¬ 
tion. We have developed several techniques for optimal extraction, including a differential method that 
eliminates low-level rogue pixels (even when no dedicated background observation was performed). The 
updated CASSIS repository now includes all the spectra ever taken by the IRS, with the exception of 
mapping observations. 

Subject headings: Astronomical databases: atlases, methods: data analysis, techniques: spectroscopic, infrared: 
general 


1. Introduction 


The Infrared Spectrograph (IRS; Houck et al. 


2004 


is one of three instruments on board the Spitzer 
Space Telescope ( Werner et al. |2004| along with two 


photometers, Infrared Array Camera (IRAC; Fazio 


et al.||2004| ) and Multiband Imaging Photometer for 


l The IRS was a collaborative venture between Cornell University 
and Ball Aerospace Corporation funded by NASA through the Jet 
Propulsion Laboratory and the Ames Research Center. 


Spitzer (MIPS; |Rieke et al.||2004| ). The IRS per¬ 
formed more than 20 000 observations over the cryo¬ 
genic mission lifetime (December 2003 - May 2009), 
corresponding to ~ 15 000 distinct targets of vari¬ 
ous kinds (Tables[l] [2}. The IRS observed between 
~ 5 and ~ 38yum in two low-resolution modules 
(R = A/AA ~ 60 - 120) and « 10 and « 37 pm in 
two high-resolution modules (R ~ 600). The main 
properties of these modules are described in Table[3] 
Most observations (« 85%) were performed in star- 
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ing mode, i.e., as single sources or groups (“clusters”) 
of individual sources. The remaining corresponds to 
spectral mappings. 


The Cornell Atlas of Spitzer/IRS Sources (CAS¬ 
SIS; http://cassis.sirtf.com), presented in 


Lebouteiller et al. ( 2011 ) (hereafter Lll), provides 


users with every low-resolution spectra observed by 
the IRS in staring mode. The pipeline performs auto¬ 
mated decisions concerning the background subtrac¬ 
tion and the choice of extraction method best adapted 
to the source spatial extent. For unresolved sources, 
the optimal extraction scales the point-spread func¬ 
tion (PSF) to the data spatial profile and provides the 
best signal-to-noise ratio (S/N) as compared to the full 
aperture extraction, with an improvement by a factor 
of two for sources < 300 mJy. Furthermore, thanks to 
the super-sampled PSF, it became possible to perform 
optimal extraction for any source position along the 
slit. The super-sampled PSF also allows users to inves¬ 
tigate complex source configurations (blended sources 
with/without extended emission component, sources 
shifted in the dispersion direction). While CASSIS 
provides the integrated spectra in such complex cases, 
the Spectroscopic Modeling Analysis and Reduction 
Tool (SMART; |Higdon et al.|2004[ Lebouteiller et al. 


20101 ) can be used for a highly-customizable manual 


extraction allowing source disentanglement. 


CASSIS represents a tool of important legacy value 
for preparing and complementing observations by fu¬ 
ture IR telescopes, in particular the James Webb Space 
Telescope (JWST). The online CASSIS database al¬ 
lows users to download spectra of publishable quality, 
and a local access to the full database is offered on 
request for large datasets. Since publication of Lll, 
several updates have been made to the low-resolution 
pipeline (see Appendix|A|). CASSIS has been used ex- 


Table 1: Spitzer/IRS observations. 



AORkeys 

Objects a 

High-res 

7192/8419 

4219/5075 

Low-res 

13565/ 16040 

10308/ 12129 

Total b 

17850/21337 

12390/ 14582 


tensively for massive dataset analysis or specific tar- 


gets (e.g., Gonzalez-Martin et al. 2015 

Brisbin et al. 

20151; Lyu et al. 2014; Sargsyan et a 

. 2014 

2012 

Brown et al. 2014; Farrah et al. 2013; Fe 

treet a . 

2013 

Gonzalez-Martm et al. 2013; Hernan-Caballero 

2012 

Hurley et al. 2012; Le Floc’h et al. 2012; Weedman 

et al. 2012). 


In the present paper, we describe optimal extrac¬ 
tion for the IRS high-resolution observations using 
a newly determined empirical super-sampled PSF. 
The two high-resolution modules use echelle spec¬ 
troscopy as opposed to long-slit spectroscopy for low- 
resolution. In the following we refer to “aperture” 
for short-high (SH) and long-high (LH) as opposed 
to “slit” for short-low (SL) and long-low (LL). The 
high-resolution modules contain 10 spectral orders 
(see Table[3] and Figure[l] for a description of the de¬ 
tectors). Staring observations work the same way as 
for low-resolution observations, i.e., a source is ob¬ 
served in two nod positions, located at about 1/3 and 
2/3 of the aperture length. With a spectral resolution 
~ 10 times higher than SL and LL, high-resolution 
observations are ideal for spectral line measurements 
and identification (and disentanglement) of narrow 
features that may be blended in the low-resolution 
spectra. For comparison, the FWHM in SH and LH 
observations, « 350 - 500 km s -1 depending on the 
spectral order, is somewhat larger than Herschel /PACS 
(60 - 320kms _1 ). Furthermore, since higher spectral 
resolution effectively results in lower signal-to-noise 
on the continuum for a given exposure time, high- 
resolution observations targeted mostly nearby bright 
sources (Table[2]). The differences between high- and 
low-resolution observations performed with the IRS 
translate into several important differences as com¬ 
pared to the pipeline for low-resolution data that was 
presented in Lll. 

We first present the pipeline steps related to the 
detector images in Sect.[2] We then explain the full 
aperture extraction in Sect|3]and optimal extraction in 
Sect.[4] In Section[5]we describe how the pipeline de¬ 
cides the best extraction method based on the source 
spatial extent. Finally, the post-processing steps at the 
spectrum level are described in Sect.[6] 


Note.— For each entry we provide the number of observations 
performed in staring mode and the total (staring and mapping). 
a Object names as given by the observer. 

^Some AORkeys were observed in both high- and low-resolution. 
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Table 2: Number of observations and distinct sources in the CASSIS atlas per scientific category. 


Category 

Low-res 

High-res 

Total 

“Cosmology” 

3/3 

0/0 

3/3 

“Cosmic infrared” 

36/36 

0/0 

36/36 

“Galaxy clusters” 

49/61 

11/12 

55/73 

“High-z galaxies” 

974/ 1258 

35/40 

993 / 1298 

“Intermediate-z galaxies” 

838 / 844 

45/46 

862/890 

“Nearby galaxies” 

462/488 

207/216 

603 / 704 

“Local Group galaxies” 

536/561 

26/27 

542/588 

“Galactic structures” 

13/13 

0/0 

13/13 

“Interacting, mergers” 

159/164 

96/97 

201/261 

“AGN, quasars, radio-galaxies” 

1394/1559 

490/538 

1549/2097 

“ULIRGS, LIRGS” 

611/632 

429 / 440 

777/ 1072 

“Starburst galaxies” 

128/133 

32/34 

153/167 

“Extragalactic jets” 

1/5 

0/0 

1/5 

“Gamma-ray bursts” 

4/4 

2/2 

4/6 

“Compact objects” 

46/59 

10/11 

49/70 

“ISM 

830/893 

279 / 290 

941/1183 

“H ii regions” 

55/59 

12/15 

56/74 

“Extragalactic stars” 

164/165 

6/6 

170/171 

“Stellar population” 

161/218 

6/6 

167/224 

“Massive stars” 

152/159 

79/81 

198 / 240 

“Evolved stars” 

1129/ 1304 

893/ 1166 

1581/2470 

“Brown dwarfs” 

250/290 

27/28 

254/318 

“Star formation” 

527/582 

95 / 102 

590/684 

“Young stellar objects” 

1532/ 1691 

336/352 

1655/2043 

“Circumstellar disks” 

2378 / 2524 

477/531 

2616/3055 

“Extra-solar planets” 

2/7 

0/0 

2/7 

“Planets” 

24/32 

16/18 

24/50 

“Satellites” 

22/23 

14/15 

24/38 

“Asteroids” 

166/170 

0/0 

166/170 

“Kuiper belt objects” 

31/31 

0/0 

31/31 

“Near-earth objects” 

12/12 

0/0 

12/12 

“Comets” 

41/54 

48/59 

77/113 

Total 

12720/ 14034 

3623/4132 

14405/ 18166 


Note.— Here we consider observations performed on sources within groups (“cluster mode observations”) as distinct observations, since several 
observations can be part of a single AOR. For this reason, the number of observations in this Table differs from what is given in Table]!] For each 
category, we provide the number of distinct sources (calculated by using a separation threshold of > 4" from other sources within the same category) 
and the total number of observations. The full list of programs with their assigned category can be found at http: //isc. astro. Cornell. edu/ 
Smart/ProgramIDs Since only one category was assigned to any given program, some sources may have a false category identification. 
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Table 3: Main properties of the Spitzer/IRS modules. 


Module 

A/AA 

Aperture size 
(") 

Pixel scale 
(") 

Order(s) 

/ lmin Amx 

(/urn) 

SL 

60-127 

3.7x57 

1.8 

1 

7.4-14.5 





2 

5.2-7.7 





bonus 

7.3 - 8.7 

LL 

60-127 

10.7x168 

5.1 

1 

19.5 - 38.0 





2 

14.0-21.3 





bonus 

19.4-21.7 

SH 

600 

4.7x11.3 

2.3 

11-20 

9.9-19.6 

LH 

600 

11.1x22.3 

4.5 

11-20 

18.7-37.2 


Note. —More information can be found at http: //cassis. sirtf.com/atlas/irs_pocketguide. pdf 


2. Detector image processing 

2.1. Individual exposures 

The CASSIS pipeline uses the Basic Calibrated 
Data (BCD) images as starting product^] along with 
the corresponding uncertainty images and the bad 
pixel mask. The BCD images are produced by the 
Spitzer Science Center BCD pipeline. We refer to LI 1 
for details on BCD images. 

Individual exposure times for SH are 6, 30, 120, or 
480 seconds. Exposure times for LH are 6, 14, 60, or 
240 seconds. There is one set of data/uncertainty/mask 
images per exposure. The mask reflects possible prob¬ 
lems identified by the pipeline. Before the exposures 
are combined (Sect. |2.2| ), the pixel masks are first com¬ 
pared over the exposures. Some pixels may have a 
lower mask value in some exposures, and therefore in¬ 
dicate more reliable values. For each pixel of the de¬ 
tector image, we select only the exposures having the 
mask values < 256. Pixels with higher mask value 
(corresponding to non-correctable saturation, missing 
data in downlink, one or no usable ramp planes, or 
pixels for which the stray-light removal or cross-talk 
correction could not be performed) are ignored in the 
other exposures. For the vast majority of cases, the 
lowest mask value is null (i.e., the pixel flux is reli¬ 
able). 


2 Tests were made with “droop” images (identical to BCD images ex¬ 
cept lacking the flat-field), but there were no visible improvements 
in the final products. 


2.2. Exposure combination 


For a given module, order, and nod position, the 
flux of each detector pixel is compared over the expo¬ 
sures. For each pixel, the exposures in which the mask 
value is relatively higher are discarded (Sect. 2.1), 
thereby selecting only the most reliable exposures 
for the combination. The final pixel flux value is the 
weighted-mean over the selected exposures. Weights 
are given according to the sequence number of the ex¬ 
posure because the first exposures are most affected 
by patterns and gradients present in the detector back¬ 
ground (Sect. 23). The relative weights are 1, 3, 5 
for the three first exposures, and 6 for the remaining 
ones. The same weights are applied when the number 
of exposures is small (e.g., weights 1 and 3 for two 
exposures). 


Since uncertainties on the individual pixels in ex¬ 
posure images may have been underestimated by the 
BCD pipeline, a simple quadratic sum of the uncer¬ 
tainties is not always accurate. For this reason, we cal¬ 
culate the uncertainties on the combined image as the 
maximum between the standard deviation of the mean 
of the pixel values over the exposures and the quadratic 
sum of the uncertainties. 


2.3. Detector background and order corrugation 

The high-resolution detectors (in particular LH) 
sometimes show light traces across or gradients. These 
artifacts manifest themselves in the exposure image in 
two different ways: 

• The background “gradient” is an unevenly dis- 
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tributed excess dark current throughout the de¬ 
tector (Fig. [2]). It affects as much as ~ 10% of the 
observations and is likely caused by a dark cur¬ 
rent residual. The detector background is gener¬ 
ally more prominent during the observation of 
bright sources, with a fraction of the source’s 
light scattered within the instrument. The detec¬ 
tor may also be affected through latency by the 
prior observation of another, bright source. The 
intensity of the background gradient decreases 
systematically with the exposure number. 

• The background “pattern”, made of diagonal 
streaks always appearing at the same locations 
(Fig. [I} arises during or after the observation 
of extremely bright sources heavily saturating 
the detector. This pattern is mostly visible in a 
handful of observations and consists of several 
clumps/traces of pixels with a non-zero level. 



Fig. 1.— LH detector image (128 x 128 px 2 ) of a heav¬ 
ily saturated source. Vertical white bands show the 10 
individual spectral orders in which the source is ob¬ 
served. For each spectral order, the wavelength axis is 
approximately vertical and the cross-dispersion axis is 
approximately horizontal. In the text we refer to rows 
and columns to describe the detector rows and columns 
in a given spectral order. The detector shows in this 
case artifacts in the form of several diagnonal streaks. 
This artifact is different from the background gradient 
presented in Fig. [2] 


The pattern appears to be related to a permanent 
bias in some pixels of the detector. 

The background pattern, when visible, is associated 
with heavily saturated sources for which most, if not 
all, pixels within the spectral orders are masked out, 
making it an irrelevant artifact to fix. In the following, 
we describe possible corrections to the background 
gradient, which is the main cause of artifacts for the 
high-resolution observations, giving rise in particular 
to spectral order tilts if not corrected. 

Removing a dedicated offset background image can 
mitigate the detector background artefact, but since the 
background gradient seems to depend on the source 
brightness, the correction is usually not satisfactory. 
Furthermore, dedicated offset background images are 
not always available or usable (Sect.[7} The difference 
between the two nod images (see Sect. |4.8| ) somewhat 
improves the removal of the background gradient but 
because of the separation in time between the two nod 
observations, a residual gradient still remains. 

In order to correct for the detector background gra¬ 
dient, one possibility is to use the dark.settlq^] al¬ 
gorithm provided by the SSC. This algorithm works 
for the LH detector only and computes a robust inter¬ 
order mean for a given row and smooths along the col¬ 
umn. It then subtracts this mean from all the data in 
the row. In this way, the inter-order region for each 
row is set to zero. The dark.settle algorithm partly 
removes the detector background but residual large- 
scale variations are often still observed. Furthermore, 
the detector background gradient is not necessarily a 
simple slope along the columns. Therefore, we de¬ 
cided to implement a custom algorithm that computes 
the smoothed 2D surface of the detector using only the 
inter-order data, and interpolates over the spectral or¬ 
ders (Fig. [2]). In practice, the spectral orders are first 
masked using a conservative mask that ensures that no 
emission from the source is accounted for. This mask 
was created specifically for this purpose. The back¬ 
ground image is then smoothed and interpolated over 
the spectral orders by means of a smoothed quintic sur¬ 
face. Note that for both dark_settle and our own al¬ 
gorithm, the calculated correction has to be performed 
on the unflatfielded image. The interpolated surface 
calculated this way is not reliable for the SH detector 
since the spectral orders in this module are too close 
to each other on the detector image, as can be seen in 


3 http://irsa.ipac.caltech.edu/data/SPITZER/docs/ 
dataanalysistools/darksettle/downloaddarksettle/ 
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Input image 



Background (model) 



Corrected image 



Background (model) 



Fig. 2.— Detector background removal for SH (£Dra, top) and LH (HD 97300, bottom). Left - Input image with the 
contrast adjusted to visualize the low-level background gradient, center - background calculated with a 2D surface 
interpolation, right - corrected image. The scale is identical for all images in a given module (set to 20% of the 
median flux within the spectral orders). The SH spectral orders are not well separated, making it difficult to determine 
the background throughout the detector. As can be seen in the top-center image, one solution is simply to extrapolate 
the background from the left side of the detector. However, considering the uncertainties on the SH background and 
the fact that the SH gradient is less problematic than the LH one, we have chosen not to correct for it (see text). 




Corrected image 


Fig. [2] Therefore we never attempt to correct for the 
detector gradient for SH. 

Despite the success of mitigating the detector back¬ 
ground gradient and significantly improving the LH 
spectra quality, our tests show that if the number of 
exposures is large enough, the combination of expo¬ 
sure images mitigates even better the gradient. The 
first exposure (in each nod observation) is indeed al¬ 
ways the most affected by the gradient, which usually 
becomes negligible after > 4 exposures. For this rea¬ 
son, we chose to apply our background removal algo¬ 
rithm to LH data only when the number of exposures 
is < 2. For a larger number of exposures, we sim¬ 
ply rely on the exposure combination with relatively 
smaller weight given to the first exposures (Sect. |2.2| ). 
The detector background is never removed for the SH 
detector, and the first exposures are simply given less 
weights. 

Despite the corrections applied above, some resid¬ 
ual emission may remain that appears as extended 
emission component in the aperture. If such a compo¬ 
nent is present, it is possible to remove it at a later stage 
when the optimal extraction is performed (Sect. [4]). 


3. Full-aperture extraction 

We describe in this section how images are used 
to perform a full aperture extraction. This extraction 
method simply co-adds the pixels in a given pseudo¬ 
rectangle (area in the detector corresponding to one 
wavelength value) to compute the flux. Since the 
flux is integrated, the presence of bad pixels anywhere 
within the pseudo-rectangle is particularly harmful. 
Therefore, bad pixels need to be identified and re¬ 
placed. The cleaning is performed on the combined 
image of all exposures since bad pixels are replaced 
using neighbors whose flux is more reliable when ex¬ 
posures have been combined. 

We use the IRS CLEAbj^] tool to substitute bad pix¬ 
els in the following order: 


• pixels with no values (NaNs) that may remain 
after exposure combination, 


pixels with a high bad pixel mask value (> 256; 
see Sect.|2.1| for the description of the corre- 


4 http://irsa.ipac.caltech.edu/data/SPITZER/docs/ 
dataanalysistools/tools/irsclean/ 
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nod 1 


nod 2 


nod 1 - nod 2 





Fig. 3.— Illustration of low-level rogue pixels being removed by differencing the two nod images. The scale is the 
same for all images. The resulting difference image can be optimally extracted if a differential PSF profile is used 
(Sect. 



sponding instrumental artifacts), 

• bad pixels and rogue pixels (i.e., pixels with a 
significant variations in their responsivity over 
time) flagged in the campaign mask, 

• pixels with a large uncertainty (> 10 times the 
median uncertainty in the image), 

• negative pixels, if > 10 times the median uncer¬ 
tainty. 

The new pixel value is calculated by the IRSCLEAN 
algorithm mainly based on neighboring pixels, al¬ 
though in some cases IRSCLEAN cannot substitute 
every eligible pixel due to clustering. Uncertainties 
and mask values are propagated for each step. The 
full-aperture extraction is performed using the stan¬ 
dard tool in SMART. Contrary to optimal extraction 
(Sect.[4]), the flux determination in the full-aperture 
extraction does not depend on the individual pixel un¬ 
certainties since the pixel values are simply summed. 
An error on the flux is ultimately calculated using the 
quadratic sum of the pixel uncertainties in the pseudo¬ 
rectangle. 

There is no background subtraction by default for 
full-aperture extraction. Although dedicated offset 
background observations may exist, they are not con¬ 
sidered in the current version of CASSIS (Sect[7]) and 
the extracted spectrum simply corresponds to the addi¬ 
tion of the source spectrum and any background emis¬ 
sion that may be present. Note that it is always possi¬ 
ble to download the full-aperture extracted spectrum of 
the dedicated background observation (if known) sep¬ 


arately and subtract it from the science source spec¬ 
trum. In this case it is preferable to subtract the back¬ 
ground at the image level to remove potentially bad 
pixels but the subtraction of the two spectra corrects 
for any emission not related to the nominal source. 

4. Optimal extraction 

Optimal extraction uses the PSF profile to compare 
to the data spatial profile in order to calculate the flux 
density (see e.g., |Home|1986| ). Optimal extraction pro¬ 
vides a spectrum with a higher signal-to-noise ratio 
when the source is unresolved. In the following, we 
describe how bad pixels are handled by the algorithm, 
how the super-sampled PSF is created, and how opti¬ 
mal extraction is performed on the data. 

4.1. Bad pixels 

For the optimal extraction, and contrary to full aper¬ 
ture extraction (Sect.[3]), the bad pixels that were iden¬ 
tified do not have to be substituted since the PSF is 
fitted to the spatial profile of the object at any wave¬ 
length. Therefore, gaps in the spatial profile are not 
problematic as long as the spatial profile is sufficiently 
sampled. In cases when the spatial profile cannot be 
reliably analyzed because too many pixels are miss¬ 
ing, the corresponding wavelength row is flagged as 
being unusable during the extraction step (Sect. |4.4| ). 
Another difference with the treatment of bad pixels be¬ 
tween full aperture and optimal extraction is that un¬ 
certainties on individual pixels are used to determine 
the flux for the optimal extraction. Therefore, bad pix¬ 
els are identified using the same steps as for full aper- 
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ture (Sect.[3]) except for the pixels with large uncertain¬ 
ties which are kept as such for optimal extraction. 

Many transient pixels, also referred to as low-level 
“rogue” pixels, are not flagged and are usually best 
corrected for by removing a background exposure. In 
the majority of observations, no dedicated background 
pointing was performed (Sect.[7]). In such cases, it 
is still possible to subtract the other nod observation 
(Fig-0 as long as one accounts for the resulting differ¬ 
ential spatial profile for extracting the flux (Sect.|4.4|). 



Fig. 4.— The super-sampled PSF for SH and LH. 
Wavelength increases from bottom to top. The black 
lines below indicate the approximate size of the aper¬ 
ture and the relative location of the super-sampled PSF 
within the aperture at the two nod positions. The 
super-sampled PSF was created on a grid with sub¬ 
pixels three times smaller than the pixel in the SH and 
LH detectors. The data beyond the secondary peak is 
that of the model PSF (see text). Such data is not used 
in standard observations in which the source lies in the 
aperture. 


4.2. Super-sampled PSF 

A super-sampled PSF, either theoretical or empiri¬ 
cal, is necessary for the optimal extraction of sources 
located anywhere in the aperture. The super-sampled 
PSF is built from mapping observations of point-like 
sources Draconis for SH and Sirius for LH) scanned 
along and across the apertures, with a step size smaller 
than the size of a pixel. We performed an iterative re¬ 
construction of the high-resolution spatial profile from 


the under-sampled data. We refer to Pinheiro da Silva 


et al. (2006) and L10 for details on the algorithm. The 


resolution on the PSF was increased by a factor of 
three for SH and LH. 


A major difference with the algorithm used for the 
low-resolution PSF described in LI 1 is that the SH and 
LH apertures are relatively small and miss a fraction of 
the desired PSF profile in any given observation. Fig¬ 
ure^ shows the relative location of the PSF within the 
aperture for the two nod positions. Since the PSF is 
never fully sampled, we built the super-sampled PSF 
“piece by piece”. For this (1) we first cut our de¬ 
sired output window (that eventually contains the fi¬ 
nal super-sampled PSF) into many overlapping sub¬ 
windows with the same size as the aperture, (2) cal¬ 
culated super-sampled PSFs for each of these output 
windows separately and (3) combined them together 
into a final PSF covering the desired range. We also 
added various steps of regularization and data filter¬ 
ing/replacement in order to improve the quality of the 
final PSF. An illustration of the super-sampled creation 
process is shown in Fig. [5] The super-sampled PSF 
provides an important improvement over the model 
created from the STINY_Tiy 2D PSF collapsed in the 
dispersion direction (green profile in Fig. [5]). . 


Our approach to the PSF is different from the 
c2d project optimal extraction (Lahuis 2007} Lahuis 


et al. 2007, 2010), the latter using an analytical cross¬ 


dispersion PSF (described as a cardinal sine function 
with a harmonic distortion component). A comparison 
of our super-sampled PSF to the c2d analytical instru¬ 
mental profile reveals a slightly narrower PSF core 
and more power in the first Airy ring (Fig. [6]). These 
differences result from the fact that the super-sampled 
PSF achieves a better resolution on the instrumental 
profile by improving on the originally low sampling of 
the PSF in the data, while the c2d PSF is calculated 
by fitting a model to the original (not super-sampled) 
data. 


A super-sampled PSF was thus created for the first 
time for SH and LH that provides the profile of a point- 
source anywhere in the aperture (Fig. [4]). While the 
sources are well-centered in the aperture in the dis¬ 
persion direction in most observations, we have also 
computed the super-sampled PSF at various positions 


5 Only the central part of the sub-windows are actually used, in order 
to mitigate edge effects related to systematic uncertainties at the edge 
of the aperture. 

t http://irsa.ipac.caltech.edu/data/SPITZER/docs/ 
dataanalysistools/tools/contributed/general 
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Fig. 5.— Illustration of the super-sampled PSF creation steps for one row of one spectral order in the SH detector 
(pixels in one row of a given spectral order have approximately the same wavelength). The thick vertical lines mark 
the region in which the PSF was calculated (the “output window”). Black segments represent input observations (i.e., 
sequence of exposures with the star located at various positions in the aperture). Blue dots show the co-addition of 
all input data (used for bad data replacement). Fractions of the PSF are calculated within “sub-windows” that are 
shown with the red dots (covering the full sub-window) and red lines (covering only the part of the sub-window that 
will contribute to the final PSF). Purple crosses indicate image data that was discarded based on our data replacement 
algorithm. The final super-sampled PSF is shown by the shifted red profile in the lower part of the plot, along with the 
STINY.TIM 2D PSF (collapsed in the dispersion direction; green), and a Gaussian fit to the core (purple). 


across the dispersion direction. Such profiles can be 
used to perform an optimal extraction of mispointed 
sources. For the CASSIS online repository, however, 
we assume the source to always be centered in the dis¬ 
persion direction. The source position along the cross¬ 
dispersion direction is a free parameter (Sect. |4.5| ). 


4.3. Spatial under-sampling 


The shortest wavelengths of SH and LH can be 
affected by spatial under-sampling. For these wave¬ 
lengths, the FWHM of the PSF is on the order of 
one pixel, which requires the use of a super-sampled 


PSF (Sect. |4.2| ) together with an accurate source po¬ 
sition. When a source shifts within a given detector 
pixel, the intra-pixel responsivity can lead to signif¬ 
icant variations in the estimated flux. Ignoring the 
under-sampling effect may result in wiggles in the ex¬ 
tracted spectrum, which is the result of the spectral 
trace (position of the source centroid in the detector 
for a given spectral order) not being perfectly orthog¬ 
onal with the detector axes. 

For the low-resolution pipeline of CASSIS, the un¬ 
dersampling was apparent for LL2 (at wavelengths 


A < 21 pm), SL2 (A < 8 pm), and SL1 (A < 14pm), by 
order of importance. The correction was performed by 
applying an empirical intra-pixel response function to 
the projected PSF on the detector grid, the corrective 
function being closer to unity with increasing wave¬ 
length. 

The smallest wavelength in SH is - lOyum, so we 
expect the under-sampling problem in SH to be equally 
important as for SL1, i.e., minor. The first release 
of the CASSIS high-resolution pipeline assumes that 
under-sampling effects in SH can be ignored. The 
smallest wavelength in LH is « 19 pm, and under¬ 
sampling problems do appear for the shortest wave¬ 
lengths (< 24pm), requiring the use of an empirical 
correction. The latter was performed the same way as 
for the low-resolution modules, with an intra-pixel re¬ 
sponsivity decreasing with distance from the pixel cen¬ 
ter and with a correction decreasing with wavelength. 

4.4. Optimal extraction kernel 

Our pipeline uses a super-sampled PSF (Sect. |4.2| ) 
and a multiple linear regression to fit the source spatial 
profile. The latter is reproduced by the combination of 
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Fig. 6.— Comparison between the CASSIS super- 


sampled PSF and the c2d analytical profile (Lahuis 


|2007|) for selected wavelengths in SH (15,20yum) and 
LH (20,30/zm). The CASSIS PSF reaches a higher 
resolution on the actual instrumental profile, resulting 
in the core being slightly narrower and the first Airy 
ring being slightly brighter. 


the super-sampled PSF itself and a large-scale emis¬ 
sion that can be parametrized (which we choose as a 
polynomial of order 0). The super-sampled PSF is first 
shifted to the source position, resampled, and finally 
scaled to the data profile. The scaling factor and large- 
scale emission are fitted simultaneously. 

As for the low resolution algorithm in L10, the 
high-resolution algorithm uses a weighted multiple 
linear regression, but within an iterative process. The 
incomplete covering of the PSF profile and the higher 
occurrence of bad pixels in the high-resolution mod¬ 
ules results in a relatively smaller number of available 
pixels as compared to the low-resolution modules. In 
the first iteration, the flux is calculated several times, 
by elevating the weight of every pixel in the core of 
the spatial profile. This gives up to 4 - 5 flux deter¬ 
minations. Outlier flux values are then flagged, and 
the uncertainty on the corresponding pixel is increased, 
since it is likely bad. In the second iteration, its weight 
is thus reduced. The iterations continue until there are 
no more outliers. In practice, only one or two iterations 
are needed. 

Several complications may occur: 

• The flux values determined using each pixel in 
a given row do not agree within errors. In that 
case, we give less weight to the pixel corre¬ 
sponding to the outlying flux value, and we per¬ 
form another iteration. 

• Pixels with no valid value are not considered for 
the fit. In order to reflect the uncertainty associ¬ 
ated with missing information, we calculate the 
fraction of flux in the valid pixels compared to 
the expected flux. The uncertainty on the final 
flux determination is then scaled up by an em¬ 
pirical coefficient inversely proportional to this 
fraction. 

• The model profile is always positive but some 
pixels can have negative values. The model is 
therefore never able to accommodate the sign in¬ 
consistency, regardless of the weights. Hence, 
rather than letting the model be biased toward 
low flux values, we replace the bad pixel with a 
null value and we increase the error bar, if nec¬ 
essary, to accommodate with the old value. 

4.5. Source finder 

The source position is first approximated from a 
Gaussian fit to the collapsed spatial profile over all 
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the wavelengths of all spectral orders. A more accu¬ 
rate position is then found through an iterative process. 
The position is varied around this first guess, and for 
each position the super-sampled PSF is calculated (i.e., 
shifted at the right position and resampled). Resid¬ 
uals are then calculated the same way as for the low- 
resolution optimal extraction ( [Lebouteiller et al. 2010 ). 
In short, the (collapsed) spatial profile of the image is 
compared to the (collapsed) spatial profile of the im¬ 
age minus the source model. The goal is to minimize 
the latter difference. The accuracy is typically less than 
1/10th of a pixel, i.e., similar to what is accomplished 
for the low-resolution algorithm. 


4.6. Extended background 


composition is available (e.g., multiple and blended 
point sources). In the website, users are encouraged 
to examine the spatial profile and validate the CASSIS 
approach. Furthermore, since in certain cases users 
may be interested in comparing the point source spec¬ 
trum to that of the extended background (or investigat¬ 
ing only the extended background), links are provided 
to download the extended background spectrum sepa¬ 
rately. Figure[7] illustrates how CASSIS disentangles 
the point-source emission from the extended physical 
background in the observation of the star HD 36917. 
The extended background is completely removed by 
CASSIS, resulting in the featureless point source spec¬ 
trum and achieving the same result as the dedicated ef¬ 
fort by Boersma et al. (2008) for this particular source. 


We refer to the extended backgound emission as the 
emission that uniformly fills the aperture and that, in 
most cases, is not associated to the nominal science 
source. The extended backgound emission within the 
aperture is either instrumental (e.g., residual from the 
detector background gradient; Sect. 2.3) and therefore 
a function of the detector row index, and/or it is phys¬ 
ical and therefore a function of wavelength. Extended 
background emission unrelated to the source can orig¬ 
inate from zodiacal emission or high galactic latitude 
cirrus clouds. In some cases, extended emission may 
be physically associated with a point source, e.g., an 
active galactic nucleus and the host galaxy. Note that 
if the science source in the aperture is found to be ex¬ 
tended, the full aperture extraction will be the default 
method and dedicated offset images need to be used to 
remove the unassociated background (Sect.[7]). 


On first approximation, the extended background in 
the aperture takes the shape of a plateau underneath the 
PSF. The multiple linear regression algorithm allows 
for a constant term to be calculated simultaneously 
with the PSF scaling factor (Sect |4.4| ). For the high- 
resolution observations however, the number of de¬ 
grees of freedom is relatively small and the spatial pro¬ 
file does not cover much of the extended background 
far from the PSF core. For this reason, a first itera¬ 
tion is performed in which the extended background 
level is derived using the multiple linear regression, 
and the derived background spectrum is then slightly 
smoothed as a function of the row index, before being 
removed from the spatial profile in a second iteration. 


The philosophy of CASSIS is to provide the spec¬ 
trum of point sources or extended sources. Apart from 
the removal of the extended background, no profile de¬ 


4.7. Uncertainties on the flux 


For all flavors of optimal extraction (Sect. |4.8| ), the 
flux uncertainty for a given wavelength element is the 
error in the PSF scaling factor. This error is domi¬ 
nated by the uncertainties on the pixels of the image 
fed to the optimal extraction core. An uncertainty is 
also calculated for the extended background, since an 
error on the latter results in an error on the point-source 
flux. Other sources of uncertainties (fringe correction, 
nod combination, calibration) are also quantified by 
the pipeline. 


4.8. Optimal extraction methods 


The optimal extraction as described in Sect. 4.4 


can 


be performed on various image products. The regu¬ 
lar method consists in scaling the PSF directly to the 
data spatial profile (Sect. |4.8.T| ) while the differential 
method uses the subtraction between two nod images 
(Sect. |4X2l ). 


4.8.1. Regular method 

Optimal extraction can be performed on the two 
nod images individually , producing two independent 
spectra which can be merged eventually into a sin¬ 
gle spectrum. The main drawback of this method is 
that few pixels are available for scaling both the super- 
sampled PSF and the large-scale emission in the aper¬ 
ture simultaneously (see Sects. |4|and |46| ). 

Optimal extraction can also be performed on the 
two nod images simultaneously. In this case, we ex¬ 
tract a single spectrum from the two images, taking 
advantage of a better sampling of the source spatial 
profile provided by the different nod positions in the 
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Full aperture extraction 




Fig. 7.— The SH spatial profile of the Herbig Ae/Be 
star HD 36917 (AORkey 11001600) is shown on top; 
it consists of a point source and a background pedestal 
emission. The lower panels show from top to bot¬ 
tom the spectra from full aperture extraction, differ¬ 
ential optimal extraction, and from the extended back¬ 
ground. CASSIS optimal extraction extracts simulta¬ 
neously the featureless point source and the extended 
background with bright emission from polycyclic aro¬ 
matic hydrocarbons at « 11.3 jam. 


detector. For this method, we assume that the PSF 
scaling factor and the background emission compo¬ 
nent will be identical for both nod images for a given 
wavelength in a given spectral order. The redundancy 
improves significantly the quality of the optimal ex¬ 
traction, in particular when bad pixels plague the im¬ 
age. The background emission is particularly better 
determined since by combining the two nod profiles, 
one can sample both sides of the PSF at the same time. 

The two methods described above are used in our 
pipeline in complementary ways, in the following or¬ 
der: 

• Optimal extraction of the two nods individu¬ 
ally to find the source position in each nod ob¬ 
servation (see explanations on source finder in 
Sect. |4.5| ). The background emission component 
is ignored for this step since it is not an impor¬ 
tant parameter for the source finder. 

• Optimal extraction of the two nods simultane¬ 
ously , using the source positions from the previ¬ 
ous step. This results in a single spectrum for the 
science object. The extended background spec¬ 
trum is saved as an optional product. 

4.8.2. Differential method 

Another method consists in subtracting the two nod 
images from each other (Fig. [5]). This subtraction is 
a reliable way of correcting low-level rogue pixels 
whose responsivity remains abnormal over a typical 
observation timescale (Fig. [3]). Since such rogue pixels 
are numerous in the detector images, a cleaning algo¬ 
rithm is rendered ineffective and in fact more harmful 
than useful. 

For the differential optimal extraction, the two nod 
images are subtracted from each other (uncertainties 
being propagated). The subtraction produces a differ¬ 
ential source profile which is scaled to the model in the 
same way as the regular method, except that the model 
is itself a differential super-sampled PSF. The latter is 
created from two normal super-sampled PSFs shifted 
in positions and inverted in flux. 

The differential method is adapted for the extraction 
of unresolved sources, even when the latter are entan¬ 
gled with extended emission. As long as the extended 
emission is physical, the difference between the two 
nod images effectively removes the background and 
produces the spectrum of the unresolved source. As 
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Fig. 8.— Optimal extraction methods for various im¬ 
age combinations. Plots are shown for a given row of 
a given order (i.e., corresponding to one wavelength 
element). Data (large diamonds) is plotted with an ar¬ 
bitrary flux density scale as a function of the pixel se¬ 
quence number along the row. In the nod 1 profile, 
one data point is missing from the profile, but the re¬ 
maining points allow for a reliable scaling. In the con¬ 
catenated nod profiles, the profiles are combined into 
a single profile and the PSF scaling factor is the same 
for both nods (although it appears different because of 
the different sampling). 


explained in Sect.[5j the differential optimal extraction 
method is not adapted for partially-extended sources. 


5. Choosing the best extraction method 


The best method for an unresolved source is un¬ 
deniably the differential method since it removes the 
background emission if present, and since it removes 
the low-level rogue pixels. If the source is not a pure 
point-like source, the differential method fails and un¬ 
derestimates the flux. 


For sources that are almost point-like source, the 
regular optimal extraction method (i.e., simultaneous 
extraction of the two nod images) provides a good al¬ 
ternative to the full aperture extraction, with a rela¬ 
tively larger signal-to-noise ratio and with a way of 
removing the background emission. Taking the ex¬ 
treme example of a source illuminating the aperture 
uniformly, the difference between the two nod images 
simply results in a image with pixel values around 0, 
with a dispersion corresponding to the S/N of the ob¬ 
servation. In such a case, the spectrum will only show 
noise. In the regular extraction of the same source, the 
PSF will be scaled to fit the flat profile as best as possi¬ 
ble, producing a spectrum with good S/N (though the 
flux calibration cannot be adequately performed, re¬ 
sulting in a wrong absolute flux and possibly a wrong 
overall spectral shape). 


Both methods of optimal extraction (regular and 
differential) are not adequate for significantly ex¬ 
tended sources because the source spatial profile can¬ 
not be modeled with a PSF. In such cases, the full- 
aperture extraction is the best method. Following the 
recent improvements of the low-resolution pipeline 
(see AppendixjA]), the best extraction method is cho¬ 
sen automatically between full aperture for extended 
sources and optimal extraction for unresolved sources 
(Sect.[5]). The source extent along the cross-dispersion 
direction is calculated by comparing the width of the 
source spatial profile to the FWHM of the PSF. Note 
that, like for low-resolution observations, we assume 
that the source is either unresolved or extended for the 
entire wavelength range. An automatic determination 
of a wavelength-dependent spatial extent (and the ap¬ 
propriate flux calibration) can only be performed in 
specific cases ( [Diaz-Santos et al.||[2010[ [2011| . The 
recommended choice of the best extraction method is 
accompanied in the website with an explanation and 
with a link to the spatial profile plot. 


When the detection level is lower than a certain em- 
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pirical threshold, the source is assumed to lie at the nod 
position and optimal extraction is chosen. The reason 
for this is that if the source is not detected, full aperture 
extraction will simply add noise while optimal extrac¬ 
tion can provide a useful upper limit on the flux, maybe 
even detecting some spectral feature. For low resolu¬ 
tion, the spatial extent determination is relatively more 
robust because the sources are small compared to the 
slit length, so the application of a PSF is more reliable. 
The spatial extent derived from low-resolution obser¬ 
vations, when available, is therefore used instead of 
the one derived from the high-resolution observation 
to determine the best extraction method. 

Examples of spectra extracted with full aperture 
and optimal extraction are shown in Figs. [9] [10| and 
El corresponding to sources of various brightness. 
As can be seen from these plots, the differential op¬ 
timal extraction generally provides the cleanest spec¬ 
trum and best overall S/N for an unresolved source be¬ 
cause low level rogue pixels are removed, but this ex¬ 
traction yields incorrect flux densities if there is any 
extension of the source. The regular optimal extrac¬ 
tion (simultaneous extraction of the two nod images) 
also provides a significant improvement in S/N over 
full-aperture extraction and is valid for marginally ex¬ 
tended sources. For extended sources, the CASSIS 
full-aperture extraction provides a significant improve¬ 
ment over the post-BCD spectrum. It should be em¬ 
phasized that optimal extraction removes the extended 
background and therefore yields a spectrum with over¬ 
all flux values smaller than the full aperture extraction. 

The choice of “best” extraction is dependent, there¬ 
fore, on real source extent beyond the spatial profile 
of the PSF. Because the measure of extent is uncer¬ 
tain, the current version of CASSIS provides simulta¬ 
neous nod extraction as the default optimal extraction 
for comparison to the full-aperture extraction. How¬ 
ever, all choices are given in CASSIS so that the user 
can compare and select a final spectrum based on the 
best estimate of source extent. 

6. Post-processing of the spectra 

6.1. Fringes removal 

Fringes originate between plane-parallel surfaces in 
the light path of the instrument. The surfaces act as 
Fabry-Perot etalons, each of which can add unique 
fringe components to the source signal. While the 
LL1 fringes are believed to be the result of a filter 
delamination discovered prior to launch, the SH and 


LH fringes originate from the detector substrates and 
probably also a filter in the IRS (Lahuis 2007} Lahuis 


et al. 2003). Fringes are removed for any extraction 


type (optimal and full aperture) with the IRS FRINGE 
tooQ The default SH and LH settings in IRSFRINGE 
are chosen to identify fringes. Fringes are looked for 
in each order individually. 


The flat-fielded (BCD) images already include a 
fringe correction, although residual fringes may re¬ 
main. Such a correction relies on the assumption that 
the fringe phase, and to a lesser extent the amplitude, 
does not vary greatly between different observations. 
We decided to remove the fringes in the unflat-fielded 
images and correct for the flat-field within the flux cal¬ 
ibration step (Sect |6.3| ). 


6.2. Spectral order overlaps 

The spectra from consecutive spectral orders in 
high-resolution observations sometimes overlap sig¬ 
nificantly. In the current version of the atlas, we do 
not attempt to combine the spectra in overlap regions 
but simply choose the order cutoffs to create a contin¬ 
uous spectrum with no overlaps. We have investigated 
a large number of sources to determine empirically the 
best cutoffs for each spectral order. The current values 
ensure that the order providing the best S/N is chosen 
for any particular wavelength range. Untrimmed spec¬ 
tra are available by request. 


6.3. Flux calibration 


The flux calibration for optimal extraction and full 
aperture extraction was performed using a relative 
spectral response function (RSRF) calculated from 76 
observations of ^Dra (including 20 before campaign 
25, .i.e., before the LH detector bias voltage was set 
to its final value). The theoretical templates and cal¬ 
ibration method is tied to the low-resolution spectral 
calibration by|Sloan et al. ( ]2014[ ) since £Dra was ob¬ 
served in both low- and high-resolution. 


Separate RSRFs were created for optimal extraction 
(regular and differential methods) and full aperture ex¬ 
traction at both nod positions. This means that each 
extraction method has been empirically calibrated us¬ 
ing a point source, so flux calibration is precise only 
for point sources. The calibration uncertainty is a sys¬ 
tematic error that depends on wavelength. 


1 http://irsa.ipac.caltech.edu/data/SPITZER/docs/ 
dataanalysistools/tools/irsfringe/ 
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7. Dedicated background offset observations 

Subtracting an offset background image can be 
a useful way to cancel out pixels with responsivity 
variations that cannot be calibrated reliably (low-level 
rogue pixels). Since the IRS high resolution modules 
were designed primarily to study emission lines, cor¬ 
rectly subtracting the underlying continuum was not 
initially considered sufficiently important to double 
recommended observing times by taking a separate 
background spectrum. Later in the Spitzer mission, 
rogue pixels developed from cosmic ray damage by 
unexpectedly strong solar flares. Only after that time 
were observers advised to include offset pointings for 
high resolution observations. As a result, a significant 
fraction of observations do not have specific offsets 
(or in some cases the offset images were not observed 
with the same exposure time as the science images). 

For these reasons and also because the identifica¬ 
tion of offset backgrounds a posteriori is not straight¬ 
forward, as there was no standard way of designing 
offset observations, we have chosen not to utilize any 
available background observations in the present ver¬ 
sion of the pipeline. For unresolved sources embed¬ 
ded in large-scale emission, both optimal extraction 
methods (Sect. |4.8| ) effectively disentangle both com¬ 
ponents. For extended sources, full-aperture extraction 
is the best method, and no background can be removed 
except a dedicated offset observation. As explained 
in Sect[3| if the observation ID of the dedicated back¬ 
ground observation is known, users can download the 
spectrum of the dedicated background observation and 
subtract it from the science source spectrum. The use 
of dedicated offset backgrounds for removal at the im¬ 
age level will be investigated in the future for CASSIS. 

8. Summary 

We present the high-resolution spectral pipeline for 
the Spitzer /IRS instrument. The corresponding atlas 
is available online at http://cassis.sirtf.com, 
complementing the existing atlas for low-resolution 
data presented in LI 1. 

High-resolution modules on the IRS are particu¬ 
larly plagued by cosmic ray hits and a particular at¬ 
tention was given to the exposure combination and im¬ 
age cleaning to remove as many bad pixels as possible. 
The pipeline produces a full-aperture extraction for 
extended sources. Unresolved sources are extracted 
with an optimal extraction using a super-sampled PSF, 


the latter created for the first time for the IRS high- 
resolution modules. Two optimal extraction methods 
are considered, (1) a method extracting the two nod 
images simultaneously and removing large-scale emis¬ 
sion that may be instrumental or physical, (2) a method 
extracting the difference of the nod images, allowing a 
complete removal of any large-scale emission and of 
some instrumental artifacts. 

We are grateful to F. Lahuis for fruitful discus¬ 
sions on the high-resolution optimal extraction tech¬ 
niques. We wish to thank again the people who con¬ 
tributed to the data reduction efforts over the IRS mis¬ 
sion. Former ISC members are especially acknowl¬ 
edged (in particular D. Devost, D. Levitan, D. Whelan, 

K. Uchida, J.D. Smith, E. Furlan, M. Devost, Y. Wu, 

L. Hao, B. Brandi, S.J.U. Higdon, P. Hall) for their 
work on the SMART software and for the development 
of reduction techniques. Moreover, our colleagues 
in Rochester (M. McCLure, C.Tayrien, I. Remming, 
D. Watson, and W. Forrest) played an important role 
in bringing additional and essential improvements to 
the data reduction used in CASSIS. This research has 
made use of the NASA/IPAC Extragalactic Database 
(NED) which is operated by the Jet Propulsion Labo¬ 
ratory, California Institute of Technology, under con¬ 
tract with the National Aeronautics and Space Ad¬ 
ministration. This research has made use of the 
SIMBAD database, operated at CDS, Strasbourg, 
France. This research was conducted with support 
from the NASA Astrophysics and Data Analysis Pro¬ 
gram (Grant NNX13AE66G) 
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Post-BCD spectrum (nod 1) 



Full aperture extraction (nod 1) 



Differential optimal extraction (nod 1 - nod 2) 




Fig. 9.— SH+LH spectrum of the bright unresolved 
star £Dra (AORkey 13349632). The best spectrum is 
provided by the differential optimal extraction. RMS 
errors are shown in dark gray and systematic errors in 
light gray. The various colors indicate different spec¬ 
tral orders. The flux density scale is the same for all 
plots to illustrate the improvement in the RMS noise. 


Post-BCD spectrum (nod 1) 
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Optimal extraction (simultaneous extraction of nod 1 and nod 2) 



Fig. 10.— SH spectrum of the relatively faint post- 
AGB star MSX SMC029 (AORkey 25646848; [Kraej 
. See Fig.|9] for the plot description. 
The emission feature at 11.3/mi and absorption fea¬ 
tures at « 13.5//m and « 16/mi are real and best seen 
in the optimal extraction versions. 


mer et alT|[2006|) 
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Post-BCD spectrum (nod 1) 





Fig. 11.— LH spectrum of the faint unresolved galaxy 
2MASX J10394598+6531034 (AORkey 33357568). 
See Fig.[9]for the plot description. Note that the source 
has no detectable emission lines. The flux scale is the 
same for all plots, illustrating how optimal extraction 
removes the background emission that is included in 
full aperture extractions. 
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Appendices 

A. Low-resolution pipeline updates 

Since publication of Lll, several updates have been made to the CASSIS low-resolution pipeline. The current 
version is labelled as “LR7”. The spectra now available through the CASSIS website (http: //cassis. sirtf.com) 
include the following improvements: 

• The best extraction method is chosen automatically between optimal extraction and tapered column extraction 
(integrating the flux within a spatial window whose width scales with wavelength), based on the source spatial 
extent. Optimal extraction is best used for unresolved sources, while tapered column extraction is adapted for 
partially-extended sources. The default spectrum shown on the main result page and the default products reflect 
the automatic choice between the extraction methods. The alternative method can still be accessed through the 
options. 

• Background subtraction for low-resolution observations was performed either by removing the detector image(s) 
corresponding to the other spectral order “by-order”) or to the other nod (“by-nod”). The presence of a contam¬ 
inating source in the nominal image and in the background image(s) is critical to constrain what background 
subtraction method is eventually used (between by-order, by-nod, or no subtraction at all). The parameters were 
adjusted so that a contamination is identified as such only when it affects significantly the source spectrum. 

• Tapered column extractions in v4 were presented with the best background subtraction based on diagnostics 
drawn from the optimal extraction algorithm (accounting for the presence of contaminating sources). If the 
source is too extended, however, it becomes impossible to disentangle the “positive” and “negative” peaks in 
the differential profile and the by-nod background subtraction is not reliable. For tapered column extractions 
of partially-extended sources, the subtraction by-order is now preferred, unless there is a contaminating source 
in the other order background, in which case no by-nod or by-order subtraction can be performed, resulting in 
what is referred to as “in situ” local background removal, which removes only the baseline to the spatial profile 
as opposed to removing a 2D background image. 

• The by-order (and to a lesser extent by-nod) subtraction sometimes resulted in a significant residual of the 
extended background emission. This mostly affects very faint sources (typically < 1 mJy) for which the source 
flux is much smaller than the difference of the background emission between two order (or two nods). The 
residual emission is now removed prior to extraction of the source profile. 

• The latest and final version of the BCD calibration is used (SI8.18.0). 

• Various improvements were made for the website, with in particular the ability to overlay the slits on archival 
images. 
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